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1 The Hydrocyclone 

The hydrocyclone is a tool used in different fields but mostly in mining 
industries to separate minerals. The minerals are put in a liquid or in a 
stable suspension of a chosen density which must be less than the density of 
the heavier minerals and greater than that of the lighter minerals [1,2,3]. 

Nowadays there is a strong interest for such tools, mostly in the processes 
where is possible to eliminate or reduce a lot the operation of grinding. This 
allows to reduce the cost of production and simplify the separation of minerals 
[4,5,6,7]. The modes of separation are two. One is the static mode, where 
the basis of the separation process is the gravitational force. The other is 
the dynamic mode, which operates by using a centrifugal force. 

The static mode can be applied when is valid the relation F g > Fp, where 
F g is the gravitational force and Fd is the resistance force to the motion of 
mineral particles into the dense fluid. The relation is strictly related to the 
mass and dimension of the particles. As a consequence, the static separation 
can be used only with dimensions of the particles greater than 4-6 mm. 
Instead the dynamic separation can be used with a dimension of the particles 
less than 0.2-0.3 mm. 

The hydrocyclone can be classified as: 

i) conic hydrocyclone (DMS); 

ii) cylindric hydrocyclone (DWP); 

iii) multi-stage hydrocyclone (TRI-FLOW). 



The conic hydrocyclone is a cylinder connected in the lower part with a 
truncated cone of an angle between 15°-20°. The discharge of materials with 
less density (overflow) is done through a tube, the so-called vortex finder, 
with center in the upper cylinder. The discharge of the heavier materials 
(underflow) is realized through a tube (apex) which is connected to the lower 
base of the reversed conical section. The principle of functioning is equal 
in all systems. The material suspended in a dense fluid (made of the heavy 
material mixed to water) is introduced as the tangent through the inlet, 
producing into the hydrocyclone a vortical movement. 

Under the effect of centrifugal forces the more heavy mineral particles 
put into the dense means will try to go to the wall of the hydrocyclone. 
In the central part along the central axis, the vortex gives the origin to 
a strong depression that allows the presence of a central nucleus (air core). 
These particles will travel along spiral trajectories and will leave the separator 
through the apex (under the effect of gravitational forces) that is in the 
lower part. Instead, the lighter particles, following the most part of the 
fluid, will tend to go in the direction of the central part, always following a 
spiral trajectory where by means of a strong ascensional kinetic field (due 
to a convenient difference of pressure applied to the two exits) will leave the 
hydrocyclone through the vortex finder put in the upper part. 

For the cylindrical hydrocyclone and for the multi-stage hydrocyclone 
there is a tangential input of the dense means to originate a vortex and a 
central depression where an air core is situated. There will be a discharge in 
the lower part. The particles that are lighter will leave by a spiral trajectory 
through the vortex finder. The heavier particles will leave after going in the 
direction of the walls through the upper discharge (sink). 

2 Confined swirling 

Because of the effect of the walls the flow in a hydrocyclone is the result of 
a movement of confined swirling (vortex) and recirculating flow, which are 
generated by the difference of the means. There are many applications of 
the movement with vortex: not only the solid-liquid separation considered 
here but also solid-gas, gas-liquid and gas-gas separation as in combustion, 
production of jet, stabilization of fluxes in a plasma, heat exchange, valves 
and turbomachines [8]. 
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The flux with vortex is the result of increasing the speed of the tangential 
component v t by using a tool called vortex generator. Vortices having enough 
force can produce in the inner part of the fluid mass some fluxes in the 
opposite direction of the vortex. The configuration of this recirculation zone 
can depend on many different factors as the distribution of velocity, the 
geometrical shape of the place where the vortex is developed, etc. 

A certain number of properties characterize the movement with confined 
vortex: 

a) the movement is 3-dimensional; 

b) the tangential component Vt of the velocity is predominant in the flux 
with the exception of the central turbulent part near the axis of the vortex; 

c) the radial component v r of the velocity is very small; 

d) the predominant component of the axial velocity v z is concentrated in a 
region close to the wall, having, close to the center, a component of inverse 
flux; 

e) the axis of symmetry can become a spiral with a dynamical behaviour in 
time; 

f) the microscopic flux is not very much related to the changes of velocity of 
the material into the system; 

g) the turbulence is always very great because the flux with vortex is in a 
system of stationary type. 

3 Properties of vortex movement 

Three are the properties that is necessary to introduce: 

1) the Rankine profile for the tangential vorticity; 

2) the vortex breakdown; 

3) processing of vortex core. 

As vortex of Rankine is intended the forced flux that can be interpreted as 
the combination of the rotation of a rigid body with contemporary presence 
of a free vortex where the action of the viscosity can be neglected. The 
rotation of the rigid body is classified as a rotational flux, instead the free 
vortex is classified as an irrotational flux. The zone where we have rotation 
of the rigid body is the central viscous region of the vortex and it is called 
nucleus of the vortex and also viscous nucleus. For a determined value of the 
axial coordinate such that the value taken by the tangential velocity is zero 
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such position is the center of the vortex which usually is different from the 
center of the place where the vortex develops [9,10,11]. 

To be able to deduce the characterized value of the Rankine vortex, one 
can follow the analysis prepared by Milne Thompson (1966). To such goal a 
plane flux is considered with a forced vortex of circular radius rj, and constant 
vorticity ( = V A v externally limited by an irrotational vortex. 

By cylindrical coordinates and avoiding the axial and azimuthal variation 
two circular trajectories are considered around the forced vortex and with 
radius T\ < r& and < r 2 indicating by and v<f, 2 the tangential velocity 
of the fluid in correspondence to the two limit trajectories. 

Applying the Stokes theorem to the flux of the vector ( = V A v through 
the circular surfaces of ray r\ and r 2 one obtains: 

j (V A v) ■ ndSi = £ vdd i = 1, 2 (1) 

where n indicates the versor of the normal to the generic surface dSi and dCi 
indicates an infinitesimal portion of the trajectory Ci of length 27rrj. By the 
development of these formulas one obtains the circulations T 1 and T 2 

= (V A v)vrl = 2vrr^ , % = 1, 2 (2) 

from which it is possible to obtain the tangential velocity for the Rankine 

vortex 

% = 2^4 (VA " )ri ' 1 = 1,2 (3) 

Remembering that the movement is realized into the plane (r, 0) the relation 
2u = V A v applied to the direction z gives 

^ = 2^r^- (5) 

Therefore the forced vortex is equivalent to a rotation of a rigid body with 
angular velocity uj. The tangential velocity v^r) represents in the cartesian 
plane a line through the origin and of angular coefficient uj with ki = u and 
k 2 = <jjr\. The tangential velocity for forced vortex becomes = = k±ri 
and for the free vortex = ^b; = ^2/^2- 
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4 Break of the vortex 



If the tangential component of the velocity, namely of the Reynolds number, 
produces a zone of centered recirculation which signs the transition from 
a condition of flux with vortex that is said to be the break of the vortex 
and depend strongly from the swirl number S. This last one is a parameter 
that indicates the measure of the tangential velocity at the input. For two 
cylinders of ray i?i and R 2 respectively that delimitate the zone of vorticity 
one has a S with the following expression 

= jg 2 pr 2 v z v r dr 
R 2 j£ prvldr • 

Sarbkaya, Escudier and Zehuder have studied the problem of break of the 
vortex. All these studies make a relation between the Reynolds number and 
the circulation number S. 

5 Precession of the nucleus of the vortex 

The recirculation zone generated by the break of the vortex is made by 
a jet with vortex in the strait direction and by the zone of inverse flux. 
These two quantities have distinct values for the axial and angular fluxes, 
sufficiently stable then exists an intermediate zone with great turbulence 
and fluctuations. This fluctuation originates the precession of the nucleus of 
forced vortex that tends to instability and moves around its symmetry axis. 
The different modes of precession of the nucleus and different stabilities are 
generally characterized using the Rayleigh number for the rotating fluid. This 
criteria says that a flux is stable if the quantity 

9 = pv^r (7) 

grows with r (forced vortex) and this quantity is stable if it is constant with 
r (free vortex) and unstable if it decreases with r. 

6 The nucleus of air 

If the energy of a movement with vortex is big enough one has the generation 
along the rotation axis of the fluid of a nucleus of air, where the transversal 
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dimension depends directly on the pressure applied to the system and the 
form of such nucleus is cylindrical only in few cases. As it was said before, 
the axis of rotation does not coincide with the geometrical axis of the hydro- 
cyclone. Namely for the nucleus of air there is no axial symmetry. Through 
the surface of the nucleus of air that is a liquid-air interphase there is a dis- 
continuity of the gradient of the radial pressure due to the presence of two 
fluids with different density. The generation of the nucleus of air is due to 
the lowering of the pressure along the center of the vortex, to simplify the 
fluid is considered non viscous as was done by Milne Thompson. Considering 
a Rankine vortex one can indicate with p\ and P2 respectively the pressure 
into the vortex region, the pressure pi with the formula 

Pi = ^ 2 r 2 + Po , (8) 
and for the potential vortex the application of Bernulli theorem gives 

P2 . 1 2 Poo / n x 

p 2 p 

where v $ = p = p^ when r goes to infinity. If the pressures are equal on 
the surface r = r b one obtains for p 

Po=Poo- p^ 2 rj , (10) 

Pi = Poo- pu 2 r 2 b (l - —5) , (11) 
1 r 2 

P2=Poo- ^P^ 2J k » ( 12 ) 

dividing by p^ the last expression and introducing the quantity 

A = pcu 2 ^- , (13) 

Too 

one obtains 

Po 2r 



( 14 ) 



Po * r b 
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The variation of the pressure in the two regions of the Rankine vortex gives 
a function of the quantity r/ro and a behaviour of parabolic type given by 
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r- {1 - A) = l A ?> (16) 

Poo * I b 



(—-1)3 = ^- (17) 



>oo rl 2 

Different formulas have been given for the radius of the nucleus of air, the 
most important is that given by Concha and Barientos as 



Ra = ~ jr- , (18) 

2fxa - Ap a 

where a and Ap depend on the geometrical parameters and \x is the viscosity 
of the fluid, a the superficial tension, a the radial component of the gradient 
of velocity evaluated in relation with R a . Ap is the difference of pressure 
through the interphase air-liquid. 



7 The dynamical model 

The model of the Hydrocyclone is 3-dimensional, when to the sink and to the 
vortex finder a difference in the pressure is applied one has the generation of 
the air core where it is stable. The equations are the Navier-Stokes for the 
conservation of the mass and of the momentum [8]. We use the finite element 
method (FEM) [12] and the k-e model [13], indicating by (u i: p,T : c,k : e) all 
the variables and taking into account that the (k, e) turbulence model the 
field, this is characterized in terms of two variables using the turbulent energy 
k which is defined as ^ 

k= -m- Uj , (19) 

e = vui ■ Uj . (20) 

In the following we will use repeated indices to indicate summation and the 
symbol a j means J^. Typical turbulent eddy velocity and length scales 

(denoted by Ui and k) can be characterized as \fk and /c L5 /e. [i t = poc^k/e 2 
is the turbulent viscosity and is directly related to the turbulent quantities 
k and e. A transport equation for k can be obtained from the Navier-Stokes 
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equations by a sequence of algebraic manipulations. This transport equation 
contains a number of unknown correlations. A second transport equation 
for e can also be derived from the Navier-Stokes equations. Application 
of a number of modeling assumptions simplifies these two equations to the 
well known equations of turbulent kinetic energy and viscous dissipation of 
the k-e model. Then the unknown variables become (ui,p, T, c, k, e) and the 
corresponding field equations are: 

u jtj = , (21) 

du ■ 

Po(-Q^ + UjUij) = ~P,i + Pofi- Po9i[Pr{T -T ) - f3 c c] + [p(u itj + u jti )j] , (22) 

dT 

Poc P (^- + UjTj) = (XT,), + p$ + H, (23) 
dc 

Po(^; + UjCj) = Poiacj), + q c + R, (24) 

+ Ujkj) = (p + + Ht® + Vt9i(—Tj + ^c :j ) - p e , (25) 

at a k a t b t 

Poift+Ujej) = (^o+— ej)j+ci^t^+ci(l-c 3 )^j(— PtTj+— (3 c c j )-poC 2 j , 

(26) 

where $ = 2e^eij^ with = ^(ui,j+Uj,i) the strain stress tensor, a is the slip 
coefficient, /3r and (3 C are the volume expansion coefficients, A is the thermal 
conductivity, gi is the gravitational force vector. After years of testing the 
k-e model, the choice of the quantity taken as empirical constants has led 
to the following recommended set of model constants: c M = 0.09, a k = 1.00, 
a e = 1.30, a T = 0.9, S t = 0.9, c x = 1.44, c 2 = 1.92, c 3 = 0.8. 



8 Turbulence Modelling 

The system of equations (21-26) is discretized by the usual finite element 
method. It is important to note that the introduction of the eddy viscos- 
ity concept makes the total effective viscosity a function of position, thus 
necessitating a stress divergence formulation. Using the approximation 

k(x,t) =<f) T K(t) , (27) 
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e(x, t) = <p T E{t) . (28) 

Within each element, the velocity, pressure and temperature field are ap- 
proximated by 

Ui (x,t) = <t> T U i {t) , (29) 

p(x,t)=^ T P(t) , (30) 

T(x,t) =6 T T(t) . (31) 

where Ui, P and T are column vectors of element nodal point unknowns and 
0, if} and 9 are column vectors of the interpolation functions. Herein the same 
basis functions are employed but cost-effective restrictions. Substitution of 
these approximations into the field equations and boundary conditions yields 
a set of equations: 

f 1 (<j>,W,U i ,P,T) = R 1 , (32) 

f 2 (<f>,U i )=R 2 , (33) 

h(<j>,6,U i ,T) = R i , (34) 

where R±, R2, and R3 are the residuals (errors) resulting from the use of the 
approximations of equation (3). The process results in the following finite 
system of nonlinear ordinary differential equations: 
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The sub-matrices M, N and Cj are given by 



M = / p^fdV , 
Jv 


(36) 


n = J v p oCp ee T dv , 


(37) 




(38) 



The remaining sub-matrices take the following form for i,j = 1,2,3. Note 
that parentheses on repeated indices indicate a relaxation of the summation 
convention: 

K M) = J v mH>jdV + J v M0,,05 + , (39) 

K l3 = [ p^dV i^j, (40) 

J V 

K l4 = f pMfdV , (41) 

J V 

^44 = J v PoCpU^dV + jf \</>j(%dV , (42) 

K m = I poUjMjdV + ( ^<f>j<^dV+ ( PoUfdV , (43) 
Jv Jv o e JV rC 

K 77 = I PoUjcfxftdV + ( ^(f>j(%dV+ I pc 2 UfdV , (44) 
Jv Jv (J £ Jv rC 

Fi = J v pfidV + j v p 9i {\ + (5T^dV , (45) 

F 4 = - [ Xe^dS + / qsOdV + / p^ddV , (46) 
Js on Jv Jv 

F 6 = I ^<pdS+ I Pt^dV, (47) 
Js <Jk on Jv 

F ?= I -7T-WS+ I c^-p^dV . (48) 
Js <r e On Jv k 

Each of the above integrals are evaluated using the isoparametric map/quadrature 
procedure. Then, using the isothermal hypothesis, the solution is obtained 
by applying the Galerkin Finite Element Method in the usual fashion. In Fig- 
ures 1-3 we show our FEM numerical results, which are in good agreement 
with the experimental data [5,6,7]. See figure captions for further details. 
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9 A simple model for the Hydrocyclone 

In this section we study a simple 3D family of maps which shows some of the 
properties of the Hydrocyclone. 

We call (p, 9, v ^-polypous the family of planar maps f p e : R 3 — > R 3 such 
that f P Vz {x,y,z) = (f^l z (x, y, z), f$ Vz {x, y, z), f$ Vz (x, y, z)), with p > 0, 
9 E [-27T, 2ir], v z E R and: 

/•(l) / \ a . ■ a x 3 y 3 

J P evA x ^y^ z ) =x pcosd + y psmd- y - — , 

(2) 

f P 0v z ( x >y> z ) = -x psin9 + y pcos9 + — - — , (49) 

fp? Vz ( x ,y, z ) = z + v z . 

In Figure 4 and 5 we show the dynamics of our map for p — 1, v z — — 1 and 
three different values of the 9 parameter. The trajectories follow a spiral and 
tend to go down in the direction of the central part. By taking a positive 
value for the velocity v z and different values for the parameter 9 we can model 
also the trajectories of particles which go in the upper part. 

Because the dynamics of the (p, 9, v ^-polypous in the z direction is triv- 
ial, we analyze in detail the (p, 9, v z )-polypous on the (x, y) plane for v z — 
(see also [14]). 

To characterize the (p, #)-polypous we can calculate the Jacobian 

D Mx, y) = ( P C0S V X \ P Sin 9 ' y l) . (50) 
Jp0y ' y; \-psm9 + x 2 pcos9-y 2 J y ' 

The trace of Df p o is given by 

Tr(Df pe (x,y)) = 2pcos9 - (x 2 + y 2 ) , (51) 

and the determinant 

det(Dfpe(x,y)) = 2x 2 y 2 - p(cos9 + sin#)(x 2 + y 2 ) + p 2 . (52) 

A direct calculus shows that the eigenvalues of the p6>-polypous are given 

by 

\i, 2 (x,y) = ^[TriDf pe (x,y))T^TriDf pe (x,y)) 2 - Adet(Df pe (x,y))] . (53) 
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The (p, #)-polypous has a fixed point in (0,0) G R 2 . This fixed point is 
hyperbolic if p ^ 1. It is asymptotically stable, i.e an attractor, if p < 1 and 
it is unstable if p > 1. The (p, 0) -polypous is differentiate but it is not a 
homeomorphism because it is not invertible: there are several anti-images of 
the origin, one of them is obviously the origin itself. 

Now we discuss in great details the very interesting case p — 1: The 
(p, 6*)-polypous has a non-hyperbolic fixed point in (0,0) G R 2 if p = 1. 
This fixed point is asymptotically stable if 9 G ] — 7r/4, 37r/4[, but 0^0. 
First we observe that when 9 = and 9 = it the eigenvalues of Df^^O, 0) are 
real numbers. For all the other values of 9 the eigenvalues are complex but 
not real and have unitary modulus. 

We call d e (x,y) = det(Df w (x,y)) and obtain 

Vde(x, y) = (Axy 2 — 2 (cos 9 + sin 9)x, Ax 2 y — 2(cos0 + sin #)?/), 

TT , , N / Ay 2 - 2(cos6 l + sintf) 8xy 
Hd e (x,y) = [ y V Q „ 2 «/ /i , • /i\ 

v ' y; y 8xy Ax 2 - 2(cos9 + sm9) 

and its eigenvalues in (0,0) are equal to — 2(cos# + sin#). The Hessian 
Hd e (0,0) is negative definite if (cos6> + sin^) > 0, thus if 9 g] - 7r/4,37r/4[. 
For these values of ^ the origin is a local maximum and by using the theorem 
of reference [14] {9 ^ 0) we have that (0, 0) is asymptotically stable. 

Consider 9 = 2n/n, n > 3. The dynamics of fw is attractive and, because 
Dj 10(0, 0) is a rotation of 9, it looks like a polypous with n branches. 

Moreover, if we add to the first two components of the (p, 9, v z )-polypous 
an oscillatory term cos (uz) we can get an elicoidal cylinder superimposed to 
the spiral dynamics (see Figure 6). So the modified (p, 9, v z )-polypous give 
rise also to the elicoidal effect, which has been recently observed experimen- 
tally [7]. 




10 Conclusions 

In this paper we have studied the turbulence and bifurcation in the motion 
of an Hydrocyclone by using a Finite Element Method based on the Navier- 
Stokes equations. The numerical results are in good agreement with the 
experimental data. We have also analyzed a simple 3D family of maps which 
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models very well the behaviour of the trajectories of the light and heavy 
particles in the Hydrocyclone. 
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Figure Captions 

Figure 1: FEM trajectories of heavy particles in the Hydrocyclone. 

Figure 2: FEM trajectories of light particles in the Hydrocyclone. 

Figure 3: FEM velocity field in the z direction for the Hydrocyclone. 

Figure 4: Spiral trajectories obtained with the (1, 9, — l)-polypous in the 
(x, y) plane. 1) 9 = 0.1 • tt, 2) 9 = 0.03 • tt, 3) 9 = 0.05 • tt. 

Figure 5: Same trajectories of Figure 4 in 3D. 

Figure 6: Elicoidal cylinder obtained with the modified (1,0.2 • n, — 1)- 
polypous with u = 0.2. 
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